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ABSTRACT 

The blazar simplified view is a new paradigm that explains well the diverse 
statistical properties of blazars observed over the entire electromagnetic spec¬ 
trum on the basis of minimal assumptions on blazars’ physical and geometrical 
properties. In this paper, the fourth in a series, we extend the predictions of this 
paradigm below the sensitivity of existing surveys and estimate the contribution 
of blazars to the X-ray and 7 -ray extragalactic backgrounds. We find that the 
integrated light from blazars can explain up to 100 % of the cosmic background 
at energies larger than ^10 GeV, and contribute ^ 40 — 70% of the 7 -ray dif¬ 
fuse radiation between 100 MeV and 10 GeV. The contribution of blazars to the 
X-ray background, between 1 and 50 keV, is approximately constant and of the 
order of 4 — 5%. On the basis of an interpolation between the estimated flux at 
X-ray and 7 -ray energies we can expect that the contribution of blazars raises 
to ^ 10% at 100 keV, and continues to increase with energy until it becomes the 
dominant component at a few MeV. Finally, we show that a strong dependence 
of the synchrotron peak frequency on luminosity, as postulated by the blazar se- 
quenee, is ruled out by the observational data as it predicts a 7 -ray background 
above a few GeV that is far in excess of the observed value. 

Key words: BL Lacertae objects: general — quasars: general — radiation 
mechanisms: non-thermal — radio continuum: galaxies — gamma-rays: galaxies 


1 INTRODUCTION 

The extragalactic sky from any direction and at all fre¬ 
quencies is filled with radiation from discrete sources and 
from a diffuse (or unresolved) component known as the 
Cosmic Background. This pervasive radiation, discovered 
only in relatively recent times, is one of the most funda¬ 
mental observables from the Universe, as it carries crucial 
information on the integrated radiation emitted over the 
entire cosmic history. The first extragalactic background 
to be detected was the Cosmic X-ray Background (CXB, 
Giacconi et al. 1962), a few years before the discovery 
of the Cosmic Microwave Background (CMB, Penzias Sz 
Wilson 1965), which is much brighter in terms of energy 
density, and is truly diffuse. The Extragalactic Gamma- 
Ray Background (EGB) was first detected by the SAS2 


* E-mail: paolo.giommi@asi.it 


satellite in the 1970’s (Eichtel et al. 1977), and for a long 
time its nature was poorly understood. The latest Eermi- 
LAT survey data are now providing strong evidence to¬ 
wards an origin mostly due to integrated radiation from 
blazars and radio galaxies (e.g. Inoue 2014; Ajello et 
al. 2015; Eornasa Sz Sanchez-Conde 2015). Blazars are 
a special type of extragalactic sources showing extreme 
observational properties, such as rapid and large ampli¬ 
tude variability, superluminal motion, and strong non- 
thermal emission across the entire electromagnetic spec¬ 
trum. Theses sources are thought to be Active Galactic 
Nuclei (AGN) that host a jet pointing almost directly 
to the observer, within which relativistic particles move 
and radiate by losing their energy in a magnetic field 
(Blandford V Rees 1978; Urry V Padovani 1995). Among 
extragalactic objects blazars are known to accelerate par¬ 
ticles to the highest observed energies and therefore are 
considered as prime candidates for multi-messenger as- 
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tronomy. Padovani & Resconi (2014), on the basis of a 
joint positional and energetic diagnostic, have suggested a 
possible association between blazars (BL Lacs) and some 
neutrino events reported by the IceCube collaboration 
(IceCube Collaboration 2014). It has also been suggested 
that blazars could be sites where ultra high energy cosmic 
rays (UHECRs) are generated (e.g. Zhang et ah 2014), 
although a firm association has so far been elusive de¬ 
spite the continuous improvement of the available data 
(P. Auger Collaboration 2014). 

Despite their low space density, compared to all 
other types of AGN, blazars are detected in high Galac¬ 
tic latitude surveys at all frequencies. In some parts of 
the electromagnetic spectrum (e.g. in the /x-wave and 7 - 
ray-ray/VHE bands) they are the dominant population 
(Giommi et al. 2009; Ackermann et al. 2011; Planck Col¬ 
laboration et al. 2011; Padovani & Giommi 2015; Arsioli 
et al. 2015, see also TevCat^). At present about 3,500 
such objects are known (Massaro et al. 2015). In the 7 - 
ray band blazars are by far the sources most frequently 
detected at high Galactic latitudes and therefore are ex¬ 
pected to be one of the main constituents of the EGB 
(e.g. Padovani et al. 1993; Giommi et al. 2006; Di Mauro 
et al. 2014; Ajello et al. 2015). 

In a series of papers (Giommi et al. 2012; Giommi, 
Padovani, & Polenta 2013; Padovani & Giommi 2015, 
hereafter Paper I, II and III) we introduced a new sce¬ 
nario for blazars based on light dilution, minimal assump¬ 
tions on the physical properties of the non-thermal jet 
emission, and unified schemes. We called this new ap¬ 
proach the blazar simplified view (BSV). By means of 
detailed Monte Carlo simulations we showed that the 
BSV scenario is consistent with the complex observa¬ 
tional properties of blazars as we know them from all 
the surveys carried out so far in the radio. X-ray, 7 -ray, 
and VHE bands, solving at the same time a number of 
long-standing issues. 

The purpose of this paper is to go beyond the sen¬ 
sitivity of current surveys and estimate the contribution 
of blazars to the CXB and EGB, building upon the sim¬ 
ulations presented in Papers I, II, and III, and applying 
the same prescriptions. 

As in previous papers we use a standard ACDM cos¬ 
mology with Hq = 70 km s“^ Mpc“^, Qrn = 0.27 and 
= 0.73 (Komatsu et al. 2011). 


2 MONTE CARLO SIMULATIONS 

In the first two papers of this series we presented the 
principles on which the BSV is built and mostly concen¬ 
trated on the statistical properties of blazars, like dis¬ 
tributions trends, average values of some important pa¬ 
rameters, such as the percentage of featureless BL Lacs, 
the observed peak frequency of the synchrotron compo¬ 
nent (z^peak)^ fhe relative fraction of ESRQs and BL Lacs, 
the amount of cosmological evolution, showing that all 
these are fully consistent with the largely different dis¬ 
tributions and values observed in real radio and X-ray 
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surveys. In Paper III we tied our simulation to the abso¬ 
lute numbers of the Fermi-2LAC catalogue and we pre¬ 
dicted the number of sources in the VHE band taking 
into account the extragalactic background light (EBL) 
absorption. In this paper, after verifying that our esti¬ 
mates are consistent with the statistical properties of the 
recently released Fermi-SLAC blazar sample, we extend 
the simulations below the limits of current surveys to cal¬ 
culate the integrated flux from the entire population of 
blazars of different types in the X-ray and 7 -ray bands. 
This requires us to explore regions of the (radio) blazar 
luminosity function not yet observed, and therefore to 
introduce some additional uncertainties arising from the 
assumptions that we must make on the slope and shape of 
the luminosity function close to the low luminosity end. 
These uncertainties are however constrained by the fact 
that here we restrict the parameters of our luminosity 
function by imposing that the predicted soft X-ray (0.3- 
3.5 keV) and hard X-ray logN-logS are consistent with 
the available data. The soft X-ray band, in particular, 
imposes tight constraints on the radio luminosity func¬ 
tion since the 5GHz radio flux density of some BL Lacs 
with X-ray flux above the limits of existing X-ray surveys 
can be as faint as ~1 mJy (Stocke et al. 1991; Rector et 
al. 2000; Massaro et al. 2015). 

Eor the reader’s convenience we briefly summarise 
here our Monte Garlo simulations (see Paper I for details) 
covering the radio through the 7 -ray bands, which are 
based on a number of ingredients, including: the (radio) 
blazar luminosity function and evolution, a distribution 
of the Lorentz factor of the electrons and of the Doppler 
factor, a synchrotron model, an accretion disk compo¬ 
nent, the host galaxy, plus a series of 7 -ray constraints 
based on observed distributions estimated using simulta¬ 
neous multi-frequency data: the distribution of Gompton 
dominance, the dependence of the 7 -ray spectral index 
on zXpgak? fhat of the 7 -ray flux on radio flux den¬ 
sity. Sources are classihed as BL Lacs, ESRQs, or radio 
galaxies based on the optical spectrum, as in real surveys. 
Readers are referred to Paper I and II for full details. 


3 BLAZARS CONTRIBUTION TO 

EXTRAGALACTIC BACKGROUNDS 

In this section we estimate the contribution to the X-ray 
and 7 -ray cosmic backgrounds by the discrete sources 
generated in our simulated survey (see Sect. 2 ). 

The contribution to an apparently diffuse extragalac¬ 
tic background resulting from the integrated emission 
from a population of discrete point sources can be cal¬ 
culated as follows: B = ^ where ^ is the 

differential logN-logS distribution, Smax is the flux of the 
brightest undetected source, and Smin is the flux pro¬ 
duced by the least luminous sources in the population at 
the largest redshift considered. 

3.1 The Cosmic X-ray Background 

The CXB has been convincingly demonstrated to be 
mostly due to the integrated light generated by accretion 
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processes onto supermassive black holes at work in AGN 
(e.g. Ueda et ah 2014, and references therein). This has 
been proven through two complementary approaches: 1 ) 
direct detection in several surveys, down to very faint sen¬ 
sitivities in the soft X-ray band (that is up to 5keV), 
and through relatively shallow surveys in the hard X- 
ray band (up to 50 keV); 2 ) AGN synthesis models, 
which estimate the contribution to the entire GXB com¬ 
bining the luminosity function and cosmological evolu¬ 
tion of the sources detected in the soft band with a pop¬ 
ulation of sources that are too heavily obscured (Gomp- 
ton thick AGN) to be detected in soft X-ray surveys (e.g. 
Ueda et al. 2014; Alexander et al. 2013). These mod¬ 
els are quite sophisticated and take into account the lat¬ 
est (radio-quiet) AGN luminosity functions, cosmological 
evolution, and parameter correlations. However, they ig¬ 
nore the contribution of blazars, which albeit relatively 
small as we will show later, is not negligible, thus intro¬ 
ducing potentially significant uncertainties and biases. 

A first attempt at estimating the contribution of 
blazars to the the X-ray and 7 -ray backgrounds was 
carried out by Giommi et al. (2006) using a semi- 
quantitative method based on the broad-band spectral 
energy distribution (SED) of blazars of different kind. 
More recently Ajello et al. (2009) made a detailed esti¬ 
mation of the hard-X-ray background based on the Swift- 
BAT survey in the 15 — 55 keV band. 

To calculate the contribution of blazars to the GXB 
up to its peak near 30 keV, we have extended our Monte 
Garlo simulations to the hard-X-ray band (30 — 50 keV) 
extrapolating the soft X-ray flux (obtained as described 
in Paper I) assuming a power law spectral shape with 
spectral index drawn from the the distribution of ob¬ 
served values in the 0.5 — 10 keV band. The latter was 
determined from a systematic analysis of all Swift-XRT 
observations of a sample of more than one hundred bright 
blazars (Giommi et al. 2015, in preparation). This was 
done both for the case where the X-ray emission is dom¬ 
inated by the synchrotron mechanism (High Frequency 
Peaked Blazars, or HBLs, with steep X-ray spectra) and 
by inverse Gompton emission (Low frequency Peaked 
Blazars, or LBLs, with flat X-ray spectra; see Padovani 
& Giommi 1995, for the original definition of LBL and 
HBL objects). 

We then calculated the logN-logS for different types 
of blazars in the soft X-ray (0.3 — 3.5 keV) band and 
compared it with observational data from the extended 
medium sensitivity survey (EMSS: Wo Iter Sz Gelotti 
2001) (see Fig. 1 ). The EMSS covers only less than 4% of 
the high Galactic latitude sky, an area that is not large 
enough to grant a sufficiently large number of FSRQs at 
fluxes larger than 10 “^^ erg cm“^ s“^ . At these high 
fluxes we estimated the logN-logS using the ROSAT all 
sky survey (RASS; Voges et al. 1999, 2000 ), which covers 
the entire sky (albeit with different sensitivities depend¬ 
ing on the exposure), and the latest version of the largest 
compilation of blazars (BZGAT5^, Massaro et al. 2015). 
We proceeded as follows: for a given X-ray flux F(0.3 —3.5 
keV) we took that part of the RASS survey where the 
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sensitivity was better or equal to F. To avoid complica¬ 
tions with the Galactic plane we only considered areas 
of the sky at |b| > 20°. We then searched for FSRQs in 
this area by cross-matching the X-ray sources with the 
BZGAT5 list. Since BZGAT5 is not a flux limited nor a 
complete catalogue, especially at faint fluxes, the number 
of FSRQs that match one of the RASS sources divided 
by the area of sky considered is a lower limit to the den¬ 
sity of sources. However, at bright X-ray fluxes (say F 
^ 10 “^^ erg cm“^ s“^ ) most of the RASS sources have 
been observed with optical telescopes and identified, and 
therefore the density estimation in this regime is reli¬ 
able. Results are summarised in Table 1 where column 
1 gives the X-ray flux F, column 2 gives the number of 
square degrees where the RASS sensitivity is better than 
F, columns 3 and 4 give the number of FSRQs and BL 
Lacs with X-ray flux larger than F that match one of 
the RASS sources within 1 arc minute in the area, while 
columns 5 and 6 give the surface density of FSRQs and 
BL Lacs. These density values are plotted in Fig. 1 as 
filled circles (black for BL Lacs and red for FSRQs^). 
The agreement between the data and our simulations is 
good. 

The simulated 30 — 50 keV logN-logS is shown in Fig. 
2. Existing surveys in this energy band have been car¬ 
ried out by Swift-BAT and by INTEGRAL (e.g. Ajello 
et al. 2009) covering a large fraction of the sky, however 
with sensitivity that is just below 10 “^^ erg cm“^ s“^ , 
i.e. orders of magnitude worse than that reached in the 
soft X-ray band. The observed surface density of Blazars 
is plotted in Fig. 2 as filled circles (green for all blazar 
types, red for FSRQs and black for BL Lacs). The agree¬ 
ment between the existing data and the simulated counts 
is good, but it is limited to the very brightest end of the 
logN-logS. Measurements at fainter fluxes will be pro¬ 
vided by NuSTAR, the first focussing orbiting telescope 
working in the hard X-ray band, which is capable of de¬ 
tecting serendipitous sources that are up to ?^100 times 
fainter than the sources detected by Swift-BAT or IN¬ 
TEGRAL. However, the field of view of NuSTAR is only 
a fraction of a degree and it will take some time before 
sufficient area of sky is covered. At the moment only one 
blazar has been reported as serendipitously discovered by 
NuSTAR (Alexander et al. 2013). 

3.2 Contribution to the 7 -ray (100 MeV - 820 
GeV) extragalactic background 

The most accurate measurement of the Extragalactic 
Gamma-Ray Background (EGB) available today has 
been recently reported by the Fermi-LAT team (Acker- 
mann et al. 2015) and covers the energy range 0.1 — 820 
GeV (see Fig. 3). 

In our calculation we must take into account the 
fact that at energies larger than ~ 50 GeV absorption 

^ The “radio-galaxies” in the figure are bona-fide blazars mis- 
classified by current classification schemes because their non- 
thermal radiation is not strong enough to dilute the host 
galaxy component even in the Ca H&K break region of the 
optical spectrum (see Paper I and II for more details). 
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Table 1. An estimate of the soft X-ray logN-logS blazar surface density from the RASS and the Roma-BZCAT catalogue (BZCAT5) 


Flux 

erg/cm2/s 

Area 
sq deg. 

No of 
FSRQs 

No of 

BL Lacs 

N(>F) 

FSRQs 

N(>F) 

BL Lacs 

l.e-12 

25600 

261 

499 

> 1.0e-2±6.3e-4 

> 1.9e-2±8.7e-4 

3.e-12 

27050 

76 

281 

2.8e-3±3.2e-4 

1.0e-2±6.3e-4 

l.e-11 

27100 

12 

91 

4.4e-4±1.3e-4 

3.4e-3±3.5e-4 



X-ray flux (0.3-3.5keV) [erg cm-^s-'] 


Figure 1. The simulated soft X-ray (0.3-3.5 keV) logN-logS 
of blazars compared to observational data. Filled circles refer 
to our estimates based on ROSAT data, while empty squares 
are from Wolter & Celotti (2001). 



X-ray flux (30-50keV) [erg cm-^S"^] 

Figure 2. The simulated hard X-ray (30-50 keV) logN-logS of 
blazars is plotted together with available observed data. The 
filled circles represents Swift-BAT measurements from Ajello 
et al. (2009), converted to the 30-50 keV band from the original 
15-55 keV band assuming a power law spectrum with spectral 
index of 1.5 for FSRQs and 2.0 for BL Lacs. 


by the Extragalactic Background Light (EBL) (Stecker 
et al. 1992) becomes important. This absorption strongly 
depends on redshift and on the 7 -ray spectral slope (F) 


in the 0.1 — 200 GeV band^. We calculate the 7 -ray back¬ 
ground at energy E as follows: 




c^-'^ebi{E,z) dN(E, S, F) 
dSdzdT 


In practice, with our simulated survey data we first build 
the differential logN-logS in narrow bins of redshift and 
spectral slopes, then we convert the integrated flux into 
flux density at energy E using the proper F, we attenuate 
the flux for EBL absorption multiplying by ^ 

and then we integrate in redshift and flux. The values of 
rebi{E,z) were taken from Dominguez et al. (2011). 

The results are plotted in Eig. 3, where the contri¬ 
bution of blazars of all types is shown as a grey band, 
representing the uncertainty due to different possible ex¬ 
trapolations of the luminosity function that are consis¬ 
tent with the observed X-ray and 7 -ray counts (see Eig. 
1 and 2). If we distinguish the contribution of the differ¬ 
ent blazars subclasses, Eig. 3 shows that BL Lacs (red 
hatched band) are the most important contributors in 
the soft X-ray band and above ~ 10 GeV, whereas ES- 
RQs (green hatched band) contribute more flux in the 
medium-hard X-ray band (~ 5-50 keV), between 100 
MeV and few GeV, and likely also at intermediate en¬ 
ergies. 

Gentral to the BSV is the assumption that there is a 
single mechanism that accelerates particles in ESRQs and 
BL Lacs and that the maximum energy that the electrons 
can reach does not depend on luminosity. In the blazar 
sequence scenario (e.g. Eossati et al. 1998; Ghisellini et 
al. 1998) instead, the maximum energy of the electrons 
is a strong function of the source luminosity. We used 
our Monte Garlo simulations to verify the implications 
of this dependence on luminosity on the contribution of 
blazars to the cosmic backgrounds. We have tested the 
correlation that defines the blazar sequence by imposing 
the following dependence of i^peak on radio luminosity: 
Log(t'leak) = axLog(!^Lj/5GHz)+ b. 

To take into account of the uncertainties in the cor¬ 
relation we allowed the parameter a to vary from - 1.2 to 
- 0.8 and set the constant b in such a way that the result¬ 
ing dependence of z^peak on is consistent with the 

famous Eig. 12 of Eossati et al. (1998). As shown in Eig. 
3 the resulting 7 -ray background is comparable to the 
observed intensity only around 100 MeV, while it rapidly 
becomes far in excess of it at energies larger then a few 
GeV. 

This excessive background value is due to the fact 
that in this scenario, HBL sources, which contribute most 


^ We have assumed a spectral break with AF = 0.5 at energies 
> 200 GeV. 
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of the 7 -ray background above ~ 10 GeV, are all con¬ 
strained in the low power end of the steep radio lumi¬ 
nosity function of blazars, on which our simulations are 
based where the source abundance is largest, thus leading 
to a very large integrated amount of 7 -ray emission. 

We note that Ajello et ah (2015) in their calcula¬ 
tion of the contribution of blazars to the EGB assume 
that the 7 -ray spectral index E (which is correlated to 
^peak) is independent of luminosity (see their eq. 12). 
With this assumption the agreement between expecta¬ 
tions and Fermi-LAT measurement is good. Had they 
assumed a strong correlation between F and luminos¬ 
ity as in the hlazar sequence^ similarly they would have 
obtained a largely different expected 7 -ray background, 
likely much larger than observed, as in the case of our 
simulations. 


3.3 The region between 50 keV and 100 MeV 

In the previous sections we estimated the contribution 
of blazars to the extragalactic backgrounds in the X-ray 
(1 — 50 keV) and in the 7 -ray (100 MeV - 820 GeV) 
bands, where data from surveys and pointed observa¬ 
tions are abundant. Galculating the contribution to the 
background at intermediate energies is subject to much 
larger uncertainties as this requires assumptions on the 
detailed spectral shape of blazars in a region of the elec¬ 
tromagnetic spectrum where only very few low-sensitivity 
observations exists. Instead of relying on unconstrained 
assumptions on the distribution of spectral shapes we 
decided to make a simple reasonable guess interpolating 
between the results in the X-ray and 7 -ray band using a 
simple geometrical shape. This is shown as a grey dotted 
band with question marks in Fig. 3. 


4 DISCUSSION AND CONCLUSION 

The blazar simplified view is a new paradigm that rests on 
a minimal set of physical and geometrical assumptions, 
and on the existence of a universal particle accelerator 
that produces the same distribution of electron peak en¬ 
ergies (see Fig . 4 of Paper I) regardless of luminosity and 
of blazar type. 

We have shown in previous papers that this scenario 
reproduces well the statistical properties of blazars in all 
the energy bands where observational data are available. 
In this paper we extended our Monte Carlo simulations to 
calculate the contribution of blazars to the extragalactic 
cosmic backgrounds in the X-ray and 7 -ray bands. Our 
results can be summarised as follows: 

( 1 ) The contribution of blazars to the 100 MeV to ~ 
10 GeV EGB is important but no larger than 50-70%. 
In this energy band other source types must contribute to 
the EGB, leaving room also for radio galaxies and star¬ 
forming galaxies that are also expected to account for a 
significant fraction of the background (Inoue et al. 2011 ; 
Di Mauro et al. 2013); 

(2) Blazars, particularly those of the HBL type, are 
responsible for 100% of the EGB above ~ 10 GeV ; 


(3) The EGB cutoff at ~ 250 GeV (Ackermann et 
al. 2015) can be fully explained by EBL absorption of 
the light that would otherwise come from distant (z > 
0.2 — 0.3) objects. This and the previous point show that 
our simulations are in good agreement with recent re¬ 
sults based on the luminosity function of Fermz-detected 
blazars (Di Mauro et al. 2014; Ajello et al. 2015) 

(4) Blazars account for about 4-5% of the GXB up to ~ 
50 keV. Assuming that the light grey band with questions 
marks in Fig. 3 is a reasonable guess, the contribution of 
blazars to the Extragalactic Gosmic Background between 
the X-ray and the 7 -ray band becomes about 10% at 100 
keV. Their contribution at larger energies increase signif¬ 
icantly and steadily until it becomes dominant around a 
few MeV (see Fig. 3). 

(5) Although blazars account for only ~ 4-5% of the 
GXB between 1 and 50 keV, this contribution is not neg¬ 
ligible and it should be properly taken into consideration 
in GXB synthesis models (see e.g. Ueda et al. 2014, and 
referenced therein). This would allow a better estimation 
of the balance between soft, hard and Gompton thick 
sources in the mix that makes up the GXB. Given that 
the contribution to the GXB increases at higher energies 
the inclusion of blazars in GXB synthesis models would 
have a significant impact on the estimation (or assump¬ 
tions) of the high-energy cutoff in Seyfert galaxies. 

( 6 ) The hlazar sequenee (Fossati et al. 1998; Ghisellini 
et al. 1998), which implies a strong dependence of the 
synchrotron peak on luminosity is ruled out by the data 
as it implies a background intensity that is at least ten 
times the observed value at energies larger than 10 GeV. 
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